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Abstract 

We continue our investigations of cavity QED with time dependent 
parameters. In this paper we discuss the situation where the state of 
the atoms leaving the cavity is reduced but the outcome is not recorded. 
In this case our knowledge is limited to an ensemble description of the 
results only. By applying the Demkov-Kunike level-crossing model, we 
show that even in this case, the filtering action of the interaction allows us 
■ to prepare a preassigned Fock state with good accuracy. The possibilities 

^— ^ ' and limitations of the method are discussed and some relations to earlier 

, work are presented. 
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; 1 Introduction 
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Cavity QED has become a standard tool of Quantum Optics. It allows the 
experimentalists to control both atomic and field variables to a high degree, 
' and consequently the system can be utilized to explore many basic features of 

matter-field interactions. In the limit when the radiation-induced time scales 
are fast enough, we can neglect the various relaxation rates affecting the system, 
and then it can be used as a laboratory realization of the well known Jaynes- 
p\ ' Cummings model. 

^ • In cavity QED physics, the model has usually been introduced with time 

independent parameters. Its utility is, however, preserved if we allow them 
to vary at a suitable rate. They should be fast compared with the relaxation 
processes of the system but slow enough not to destroy the mode structure of 
the cavity. In our works [Q and 0, we have considered the new possibilities 
offered by such a model. 

Because the J ay nes- Cummings model separates into decoupled two-level sys- 
tems, we advocated in reference pQ that it may be used to design a filter mech- 
anism to shape the photon distribution in the cavity into a desired form. There 
are earlier investigations into the effect on the cavity field caused by recording 
of the atomic state, see for example [H] and 0], but we explore the possibilities 
deriving from time dependence of the coefficients. In reference [2] , we consider 
the possibility that the desired effects may be derived simply from the shape of 
the eigenmodes of the cavity field. 
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In our earlier works, we have assumed that we record the state of every 
atom exiting the cavity after an interaction. In practical situations, however, 
the environment rapidly destroys the coherence between the atomic states af- 
ter leaving the cavity, even when no measurement is performed. In this case 
the state of the field is reduced to one of its possible configurations without us 
knowing which one has occurred. We have a situation of nonselective measure- 
ments. After the passage of m atoms through the cavity, our uncertainty about 
the actual state of the cavity becomes large. This, however, is just the situation 
to be described by a statistical ensembles of fields, each member occurring with 
a probability determined from the sequence of atomic reduction events which 
brought it about. In this paper we consider the evolution of the field in such a 
situation. 

The approach may be useful even if we, for some reason, do not want to 
record the results of the individual interaction sequence. In particular, we ad- 
dress the question of creating a pure Fock state of the field in the cavity. This 
problem has been discussed before but all methods tend to become difficult 
if one aims at large quantum numbers. Our results may be considered to be a 
contribution to the discussion aiming at the creation of Fock states. 

In sectional we set up the problem and introduce the basic concepts. In 
section |3 we specialize the Hamiltonian to the solvable Demkov-Kunike model 
and discuss its properties in its adiabatic and non-adiabatic limits. In the former 
case we show that it is, ideally, possible to reach arbitrary n-states of the field. 
Special cases of the latter situation have been discussed earlier 0] and and 
hence we only make a numerical investigation of the properties of the result 
specifically addressing the question of obtaining a Fock-state. Here the concept 
of trapping states becomes central, and we suggest a way to improve the 
method by velocity selection of the atoms sent into the cavity. In section 0] we 
conclude the presentation with a discussion. 



2 Adiabatic Jaynes-Cummings model 

We consider a system described by the Jaynes-Cummings model |HJ, which is the 
simplest formulation of the interaction between matter and a quantized field. 
The model assumes that the time of interaction is much shorter than other 
time-scales in the system, such as atomic and cavity field relaxation times. This 
means that atomic and field losses can be neglected during the interaction. The 
Hamiltonian for the full system is (h = 1) 

Hjc = fia t a+^er 3 + g(a}a~ + acr+) , (1) 

where the ct:s are the ordinary Pauli matrixes and {a, a^} are the Boson opera- 
tors of the cavity mode. Within this rotating wave approximation, the number 
of excitations is conserved; the operator N — a^a + \oz is a constant of mo- 
tion. We therefore define the interaction Hamiltonian Hi = Hjc — CIN, which 
separates the dynamics into two-level manifolds within the bare states. 
The state of the whole atom-field system can be expressed as 

oo 

I*) = c o a_(0)|0, -) + J2°n [a+(n)\n - 1, +> + o_(n)|n, ->] (2) 

n=l 
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and after the passage of the atom, the state of the system is determined from 
a± (n) . The coefficients a± (n) are given by the Schrodinger equation 



dt 



a+(n) 
a_ (n) 



Auj r - 

r- AbJ 



a+(n) 
a_(n) 



(3) 



Here the detuning Aw = cj — O and the coupling 5 are allowed to depend on 
time. Initially the atom is taken to be either in its upper or lower level 



Case (a) 
Case (b) 



Kl = i, 

a° + = 0, 



a a _ = 
l«° I = 1 



(4) 



and the initial field distribution is determined from c n . 

After the interaction, the atomic state is measured, and from equation (J2j 
it follows that the initial state of the field is modified by a+(n) or a£?(n), 
depending on whether the atom is found in its upper |+) or lower |— ) state. 
After 771 atoms have passed through the cavity and been measured, the field 
is multiplied by a series of a± (n):s. It is important to note, that the number 
of atoms detected in their upper and lower levels is not enough to tell what 
the state of the field is. One also needs the order of upper and lower detection 
events. This follows from the fact that, if the initial and final atomic states 
are different, then not only will the field be modified by a±(n), but the photon 
distribution will also be shifted by one unit. For example, if an initial lower 
level atom is found in its upper level after the interaction, the distribution of 
the field is Px t +(n) — N\a+(n+ l)| 2 P (n+ 1), where Po(n) is the initial photon 
distribution and N a normalization constant. 

Let k be a vector denoting the sequence of measured atoms defined by: 
kj = ±1 if the j:th atom's initial and final measured states have been flipped 
according to |=p) — > |±), and kj = if they are the same. Then k uniquely 
determines the final photon distribution, provided that the initial state of the 
field is known. 

If we consider Case (b) in equation and have, for example, m = 3 and 
k = (0, 1, 1), that is the first atom is found in |— ) while the other two atoms are 
found in | +) , then the final photon distribution will have the shape 



P m ,k(n) = |«+ (n + 1)|>+ (n + 2)\ 2 \a™(n + 2)\'P (n + 2). 



(5) 



If v = J2j kj , then the photon distribution after m atoms can be written, in 
a compact form, as 



(6) 



where A m ^(n) is the appropriate sequence of "filter functions" \a± (n)\ 2 and 
TVk is the normalization constant. If the field was initially in a Fock state, then, 
whatever the filter functions are, the final field will also be in a Fock state, 
shifted by an amount —v. 
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If the initial conditions are as Case (a), the probability for measuring the 
atomic states |±), given a normalized photon distribution P = P{n), is 

P(+\P) =E„K(n + l)\ 2 P(n) 

(7) 

P(-\P) = l-P(+\P). 

By repeatedly applying the equation above, we generate the probability P(k|Po) 
of detecting the sequence k, if the initial photon distribution is Pq (n) . Note that 
this probability equals the normalization constant in equation © 

P(k|P ) = N k . (8) 

Now suppose that we perform a nonselective measurement of the atomic 
state, i.e. that is we do not record the final state of the atom. Then the field is 
described by a density operator, obtained by tracing over the atomic degrees of 
freedom 

p field = Tlatom{p) =J2cnC^ [ a ^(n) a ^*(n')|n-l)(n'-l|+a!°(n)a! *(n')|n)(n'| 

n,n' 

(9) 

Note that, if we only consider the diagonal elements, P(n) = Pnn ld , we do not 
need the phases of the a±(n):s. This P(n) is the ensemble photon distribution 
as will be explained below. For Case (a) in equation (0J, the initial photon 
distribution is Po(n) = |c„+i| 2 and, after m nonselective measurements, the 
ensemble photon distribution is obtained from the recurrence relation 

P m (n) = \a™(n + l)\ 2 P m ^(n) + |o~ '(n)\ 2 P m _i(n - 1). (10) 

A special case of this equation has been studied before, see 0] and It is 
of essential importance to point out, that in general P m (n) ^ P m ,k("-), for an 
arbitrary k. The distribution given in l|10f) is the ensemble average over all 
possible outcomes P m ,k(n) 

Pm(n) = ^(k|Po(n))P m , k (ra) = ^ A mM {n)P (n + v) = P m , k (n), (11) 

Allk Allk 



where we have used equations © and Q in the second step and P m ,k(^) defines 
the ensemble average over all possible k:s. If we consider cases with a fixed series 
of outcomes defined by a given k, we have 

P m ,k'{n) =P m , k (n), Vk, k' P m (n) = P m , k (n). (12) 

This happens if for example P(k|Po) = 1 for any k; then this outcome has unit 
probability and is the only possible one, all other outcomes have zero probability. 
The only possibility that P m (n) will describe a Fock state exactly, is if all 
Pm,k(?i) represent the same Fock state. In other words, if P(n) describes a Fock 
state, this agrees with the actual state of the field. 

With the above formalism, expectation values (■■■) m with respect to the 
distribution P m {n) are related to expectation values (...) m ,k for the distribution 
Pm,k(^) according to 

(."•> m = <^W- (13) 

~ 2 

For example, the variance An m of the distribution P m (n) is given by 

An m = (n 2 ) m , k - ((n) m , k ) • (14) 
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3 The Demkov-Kunike model 

The Hamiltonian in equation J3J is in the most general form of a 2 x 2 hermi- 
tian matrix [9|, and has been widely investigated. For a review of analytically 
solvable time-dependent two-level systems see Here we only consider the 
Demkov-Kunike model, in which the two parameters of the Hamiltonian are 
given by 

^ = £ + £ tanh(^) 

(15) 

g(t) = 5o sech(£) . 

If we consider Case (a) of equation 0} and integrate over time (— oo, +oo), the 
filter functions become 



K(n)\ 



2 



_ sinh irT(E + y/E%-g%n) sinh ttT (g -^/ffg -ggn) - 

cosh kT(E+E q ) cosh ■kT^E-Eo) > 



cosh 7rT (e+^/eJ — SqwJ cosh ttT^E— \J E 2 — g 2 nj 
cosh irT^E+Eo) cosh 7rT(B-Bo) 



|a_»(n)| 2 = 1 - \a¥(n)\< 



(16) 

These expressions contain the three dimensionless parameters 

Ax = ET, A 2 = E T and rj = g T. (17) 

Normally the adiabatic limit arises from large T, or equivalently that, at least 
one of E or Eg becomes large. In the non-adiabatic limit E w Eq s» 0. In the 
following we consider these two regimes separately. 

3.1 The adiabatic limit 

With the parameters introduced in (II 711 the filter functions may be written as 

cosh7r2A 2 , cos^rgyg^Af) 



|a3°(n)r - 



i cosli 7T2A 2 I V v ' ±J A > A 

2cosh-n-(Ai+A 2 ) cosh7r(Ai-A 2 ) 2 cosh 7r(Ai +A 2 ) cosh 7r(Ai -A 2 ) 1 2 

i oa cos ( 7r2\/77 2 n— A? ) 
cosh7r2Aj | V V ' 2 I A A 

2cosh7r(Ai+A 2 )cosh7r(Ai-A 2 ) 2 cosh tt(Ai +A 2 ) cosh tt(Ai - A 2 ) - 



For no photons, n — 0, we have |aS°(0)| = 1. For 



(18) 



n> {i)> (19) 

the last term is an oscillating cosine-term in n. In the adiabatic regime, either Ai 
or A2 becomes large, and since the denominator grows exponentially in these 
parameters, the oscillating cosine-term can be neglected. This means that, 
for photon numbers n > (A2/77) 2 , the filter functions will become constant, 
determined by the first terms/term in . Then, if the parameters are chosen 
such that A 2 /f] « 1 and at least one of Ai or A 2 is large, we have to a good 
approximation 

1; n = 

\^:i»f = { , (20) 
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where < k < 1. The number n depends on the parameters Ai and A2. For 
example, if Ai = ±A2 then k = 0.5 or if Ai < A2 (the levels cross) we have 
K = 0, while A2 < Ai (no level crossing) gives k = 1. If there are non-adiabatic 
contributions, it is possible to have a k that differs from the three cases given 
above. An important point is that, as long as the process is adiabatic, the 
K = 0, 1 cases are not sensitive to the parameters of the problem. On the other 
hand, the n = 0.5 case is sensitive to the condition |Ai| = | A_2 1 - The filter 
function l|20[l with k = 0.5 can also be achieved by using a constant detuning 
Au and an asymmetric field coupling g(t) in the Hamiltonian see 
With the result l|2()|l. the equation (jlOJI becomes 

P m {n) = /tP TO _i(n) + (1 - K)P m -i{n - 1). (21) 

This equation can be solved by introducing a generating function 

00 

G m (z) = J2 znp m(n). (22) 

The equation (|21|l then becomes 

G m (z) = [K+(l- K )z]G m . 1 (z), (23) 

having the solution 

G m (z)=[ K +(l- K )z] m G Q (z). (24) 

With the field initially in a vacuum, Po(n) — <5„.o, we find after m atoms the 
state 

G m (z) = [ K +(1- k)z]" 1 = ]T ( ™ J n m ~ n {\ - n) n z n . (25) 

n=0 ^ ' 

As the coefficient of z n is the photon probability distribution, we find 

P m (n) = ( ™ "(l-«) n . (26) 

Pm{n) is the form of a Binomial distribution with a variance 

An m = tor(1 - k). (27) 

If k = we note that we have P m {n) = 5 r ^ m . From the discussion above, 
we know that if the ensemble state of the field is in a Fock state then the real 
state of the cavity field is in the same Fock state. In other words, after m atoms 
has passed through the cavity, the state of the field, initially in vacuum, will be 
|m). A similar idea for producing a Fock state has been presented in |12j . 

Physically, the process can be understood from the fact that as each atom 
traverses the cavity, it is tuned across the field-mode frequency in an adiabatic 
manner such that it exits the cavity in its lower state. Then, after the passage, 
the mode will contain one extra photon. On the other hand, if we consider the 
initial condition according to Case (b) in equation instead, the mode will be 
left with one photon less. This holds for all initial distributions of the field 

P m (n) = Po(n — m) for Case (a) 

(28) 

P m (n) = Po(n + m) for Case (b). 
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Obviously, Case (a) will heat up the cavity mode, while Case (b) cools it down. 
If lower level atoms are sent through the cavity, the field mode will eventually 
end up in the vacuum. The process can then be used as a "field eraser" |13| . 
Note that the result in equation 1(28(1 agrees with the one in 

3.2 The non-adiabatic limit 

If the atomic and field-mode frequencies are on resonance (E = Eq = 0), the 
filter functions are given by 

\a+ (n)\ 2 — cos 2 (irrjy/n) 

(29) 

\a°°(n)\ 2 =sin 2 {tttj^H) , 

where the initial condition has been chosen according to Case (a). The ensemble 
distribution (|10|l . with the filter functions as in equation 1(29(1 . has been inves- 
tigated in 0] and so we will not study it in great detail here. However, we 
specifically investigate its prospects to speed up the production of a Fock state. 
Figure ^ shows how the distribution P m (n) , initially in a coherent state with 
n = 47, evolve as the number of atoms m increases. The filter functions are 
given by 1)29(1 with r\ — 1 . We see how the distribution builds up around photon 
numbers n' = 35, 48 and 63. These n':s fulfill the trapping condition [7] 

y/(n' + l)r, = q, 9 = 1,2,3,.... (30) 

The corresponding states { |rj.') }, fulfilling the trapping condition (|30|) . are called 
trapping states and are labelled by (n',q). If the field is in a trapping state, 
then, since |aSf(n' + 1)| 2 = 1, all atoms will exit the cavity in their upper level 
and the state of the field is just multiplied by a phase-factor. The integer q 
is the number of the corresponding maximum of |a^°(n')| 2 , in other words the 
index of the trapping state. For a given trapping state (n',q), there are q — 1 
other trapping states with photon numbers smaller than n' . q also gives the 
number of Rabi cycles during the interaction |14j . When the atomic velocity 
can be controlled and consequently the effective interaction time T oc 77, we may 
change the number of Rabi cycles q in the trapping state (n', q), since T ~ q. 

Since \a°°(n' + 1)| 2 = for the trapping states, it is clear from (jTU|> that the 
states \n') are not coupled to the states |n'+l) . This results in a separation of the 
space {\n)} into disconnected blocks; population within one block will not leak 
into another. If q is small, the distance between two consecutive trapping states 
is large, and may exceed the width of the initial photon distribution. Then, if 
this distribution lies in one such block between two maxima of |a^°(n')| 2 , it will 
eventually end up in the trapping state with the largest photon number n' . In 
this way, the situation allows us to create any Fock state. However, using this 
method, the preparation of Fock states with large photon numbers is difficult 
for two reasons: 
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1. When the photon distribution "approaches" the trapping state, 
more and more atoms will exit the cavity in their upper state 
and leave the field approximately unchanged. This follows from 
the small curvature of the of the filter functions (|29() around one 
maximum with a large n' and a small q. The larger the photon 
number n' the smaller the curvature, and the number of atoms 
needed for creating the Fock state is growing rapidly. This means 
that relaxation of the field may not be negligible. 

2. The error in the atomic velocity, or in the interaction time T, prop- 
agates to n'. When this error is of the order of unity, the trapping 
condition is clearly destroyed. This undesired effect manifests it- 
self for large photon numbers, since the uncertainty in n' grows 
linearly with the photon number. 

The small curvature of the filter functions slows down the process, so if 
the curvature could be increased as the distribution approaches the trapping 
state, the number of atoms needed would be reduced. This can be achieved by 
changing the atomic velocity, and consequently the trapping state, in such a 
way that (n',q) — » (n',q r ), where q' > q. 

In figure we have plotted P TO (10), when the initial state was a coherent 
one with n = 4. The solid line (a) gives the probability when all atoms have the 
same velocity, such that n' = 10 fulfills the trapping condition (|30(1 with 5=1. 
For the solid line (b) , the velocities of consecutive atoms are chosen such that 
q is increased by one unit for each atom. Starting with q = 1 it means that the 
first atom has a velocity in which the state (n',q) = (10,1) fulfills the trapping 
condition 1)30(1 . and the velocity of the second atom is changed such that (10, 2) 
instead satisfies the condition. The trapping states thus occur according to 

(10,l)->(10,2)->(10,3)->... (31) 

for each atom. This sequence may not be the most efficient one, but it clearly 
shows how the process is speeded up. The other two lines (dot-dash) in figure |2] 
describe the same processes as the ones above, but with a 2 % Gaussian random 
error in the effective interaction time T. The error in T clearly "washes out" 
the trapping condition as mentioned above. 

4 Conclusion 

In this paper we have continued our investigation into the possibilities offered by 
letting the parameters of an atom-cavity system depend on time. It is essential 
that this variation is slow enough to retain the identity of the cavity modes 
and atomic states; adiabatic changes are assumed. On the other hand, the 
interaction times considered have to be short enough that relaxation processes 
can be neglected. 

In our earlier publications we have investigated how such time dependence 
combined with observation of the state of the outgoing atoms may be used 
as a photon state filter. Choosing a suitable time variation, desired photon 
states may be achieved. In this paper we continue the investigation by assuming 
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reduction of the atomic state but no recording of the result. This may take place 
spontaneously or by special design. The result is, however, that knowledge about 
the cavity state is assembled into a density matrix, representing the ensemble of 
events compatible with the events having occurred. In general this conveys less 
information about the state than a selective measurement does. However, if the 
method can be used to prepare a pure Fock state, the result becomes uniquely 
defined. 

The analysis utilizes a Jaynes-Cummings model with time dependent coef- 
ficients. This splits the problem into blocks that can be analysed separately. 
Using the well known level crossing model presented by Demkov and Kunike, 
we may look separately at the adiabatic and non-adiabatic limits of this model. 
Our numerical investigations suggests some scaling of the process: 

In the adiabatic limit it takes n atoms to prepare the n :th Fock state. In the 
non-adiabatic limit and with the trapping state fixed, the number of atoms grows 
with n 2 which may conflict with the requirement of no relaxations occurring. 
If, however, the consecutive atoms utilize consecutive trapping states, i.e. their 
velocities are fine tuned before entering the cavity, the number of atoms may 
be reduced to n again. 

We have not been able to solve the general problem: given that we want 
to reach an arbitrarily chosen final photon state, is there a time dependence 
which guarantees this? We have only illustrated the utility of the method for 
the preparation of a selected Fock state. 
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Figure 1: The probability distribution P m (n) is plotted versus the the number 
of atoms to. Here the dimensionless parameter r\ = 1. The distribution is built 
up around the photon numbers n = 35, 48 and 63. These photon numbers all 
fulfill the trapping condition 1)30(1 . 



Figure 2: This figure shows how the probability P m (10), i.e. having n = 10, 
evolves with the number of atoms to. For the solid line (a), the velocities of 
the atoms are all the same and chosen such that n = 10 fulfills the trapping 
condition l|30|) with q — 1. For line (b), the velocities are changed for each atom 
in such a way, that the trapping states (n,q) occur in the sequence (10, 1) — > 
(10,2) — > (10,3)..., for consecutive atoms. The two dot-dash lines are the same 
as the solid lines, but with a 2 % Gaussian random error in the atomic velocity 
and consequently the interaction time. The spread in velocities clearly destroys 
the trapping effect. 
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